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ABSTRACT 



Determining the properties of remote globular clusters and elliptical galaxies using 
^ \ evolutionary population synthesis requires a library of reliable model stellar fluxes. Em- 

■^ I pirical libraries are limited to spectra of stars in the solar neighborhood, with nearly 

solar abundances and abundance ratios. We report here a first step towards providing 
a flux library that includes nonsolar abundances, based on calculations from first prin- 
ce ■ ciples that are calibrated empirically. Because the mid-ultraviolet spectrum of an old 
JtL . stellar system is dominated by the contribution from its main-sequence turnoff stars, 
^^ , we have started by modeling these. We have calculated mid-ultraviolet spectra for the 
^«^ \ Sun and nine nearby, near-main-sequence stars spanning metallicities from less than 
^^ \ 1/100 solar to greater than solar, encompassing a range of light-element abundance 
enhancements. 

We first determined temperatures of eight of the stars by analyzing optical echelle 



o 
o 

Q^' spectra together with the mid-ultraviolet. Both could be matched at the same time only 

Q . when models with no convective overshoot were adopted, and only when an approxi- 

^ , mate chromosphere was incorporated near the surface of relatively metal-rich models, 

j^ i Extensive modifications to mid-UV line parameters were also required, notably the man- 



ual assignment of approximate identifications for mid-UV lines missing from laboratory 
linelists. Without recourse to additional missing opacity, these measures suffice to re- 



Vh ' produce in detail almost the entire mid-UV spectrum of solar-temperature stars up to 



one-tenth solar metallicity, and the region from 2900 A to 3100 A throughout the entire 
metallicity range. Ramifications for abundance determinations in individual metal-poor 
stars and for age-metallicity determinations of old stellar systems are briefiy discussed, 
with emphasis on the predictive power of the calculations. 
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1. Introduction 

To characterize the metallicity and age of an old stellar system, iron abundance [Fe/H] (the 
logarithm of the iron-to-hydrogen abundance ratio with respect to its solar value) and the stellar 
effective temperature T^s must be determined together for its main-sequence turnoff (MSTO) stars. 
These are dwarfs and subgiants of types F and early G with 5750 < T^s < 7000 K. Their spectra 
are dominated by Fei lines, whose strength increases as T^s declines. For both metallicity and age, 
then, a reliable derivation of Tpg is necessary. 

This can be problematical even for an individual star. The photometric colors often used to 
derive Tefj are sensitive to reddening by the interstellar medium and to the modeling of convection 
(Castelli, Gratton, & Kurucz 1997), as is the more recently developed infrared flux method (Alonso, 
Arribas, & Martinez- Roger 1996). Methods using Balmer-line profiles are reddening-independent 
but still very model-dependent (Magain 1984; Fuhrmann, Axer, & Gehren 1994). Methods relying 
on the excitation equilibrium of Fe i are reddening-independent and model-independent if lines are 
weak, but are susceptible to systematic errors in laboratory measurements of gf-values (Blackwell 
et al. 1982) and to possible departures from local thermodynamic equilibrium (LTE; Thevenin & 
Idiart 1999). 

Currently, Tgfj values derived for individual MSTO stars show uncomfortably large spreads: 
values of Carney et al. (1994) and Fulbright (2000) average 100 K - 200 K lower than those found by 
King (1993) and Gratton, Carretta, & Castelli (1996), with those of Alonso, Arribas, & Martinez- 
Roger (1996) and Blackwell & Lynas-Gray (1998) in between. Not only is overall metallicity 
affected, but also the ratios of the abundance with respect to iron of certain elements. Among 
these are the [Li/Fe] and [0/Fe] ratios so critical to uncovering the nucleosynthesis history of 
the early universe (Spite et al. 1996) and deducing the age of globular clusters from the MSTO 
luminosity (VandenBerg & Bell 2000). When Tgff alone is raised by 100 K, modeled values of [Li/Fe] 
change by -h0.09dex (Spite et al. 1996); [0/Fe] changes by -0.09 dex, -F0.03dex, or -h0.16dex 
depending on whether the high-excitation Oi 7775 A triplet, the [Ol] 6300 A ground-state line, or 
OH molecular lines near 3100 A are used (Kraft 2000). Not surprisingly, the oxygen abundance in 
extremely metal-poor stars remains uncertain, as summarized during a workshop devoted to the 
problem (Kraft 2000; Balachandran 2000; Peterson 2000; Lambert 2000). King (1993) and Kraft 
(2000), among others, have suggested that a higher Tefj scale for metal-poor turnoff stars might 
eliminate the systematic differences in [0/Fe] seen from its different diagnostics. 

For all these reasons, it seems important to reduce uncertainties in the determination of T^fi 
in MSTO stars. Fitting the slope of their UV flux distribution might help, since the fluxes of F-G 
stars peak in the optical so that UV fluxes vary rapidly with Tcs . Although fitting a UV slope is 
also reddening-dependent, reddening should be very low for the brightest metal-poor MSTO stars, 
which are typically within 50 pc (e.g., Reid 1998). 

Reliable calculations of MSTO mid-UV spectra would be especially valuable in constraining 
the age and metallicity of a spatially unresolved old stellar system such as an elliptical galaxy or 



extragalactic globular cluster. For the mid-UV flux of an old stellar population is dominated by 
turnofF stars (Dorman, Rood, & O'Connell 1993), but the light redward of 4500 A is dominated by 
cool giants in a 17Gyr population of solar metallicity (Worthey 1994). To date, however, mid-UV 
studies of galaxies and globulars (e.g., Spinrad et al. 1997; Ponder et al. 1998) have largely been 
empirical, comparing extragalactic spectra against observed spectra of individual F and G stars 
(Fanelli et al. 1990). Thus the metallicity of the extragalactic systems has generally been assumed 
to be solar, although age is very sensitive to this assumption (Heap et al. 1998). Furthermore, 
nonsolar abundance ratios have not been included, despite the fact that abundance enhancements 
of light elements such as magnesium are generally seen in metal-rich extragalactic systems (Peterson 
1976; O'Connell 1976; Worthey, Faber, & Gonzalez 1992; Henry & Worthey 1999). 

Modeling the mid-UV region is difficult, though, because of the dramatic increase in line 
absorption at short wavelengths. Line blending in the solar spectrum suppresses the optical con- 
tinuum blueward of ~ 4500 A (Kurucz et al. 1984), and becomes even more severe in the near- 
and mid-UV. Especially in the UV, many lines are "missing" - although they appear as absorption 
features in spectra of solar-type stars, they lack laboratory identification, and are thus not matched 
by spectral synthesis calculations (Kurucz & Avrett 1981). As a result, the normalization of the 
solar UV continuum is controversial (Balachandran &: Bell 1998; Boesgaard et al. 1999; Israelian 
et al. 2001). 

With few exceptions (e.g.. Heap et al. 1998), calculations of the mid-UV region are based on 
Kurucz model flux distributions. Allende Prieto & Lambert (2000) compared these to International 
Ultraviolet Explorer (lUE) mid-UV spectra of nearby metal-poor halo stars, illustrating that T^s is 
generally recovered when parallaxes are known and reddening negligible. Lotz, Ferguson, &: Bohlin 
(2000) have used Kurucz fluxes as generated by Lejeune, Cuisinier, & Buser (1997) to examine the 
effect of abundance on the mid-UV light of old stellar populations, but found that the UV fluxes 
of cool stars were poorly reproduced. 

This is not surprising given the known limitations of the opacity distribution functions on 
which the Kurucz flux distributions are based. They statistically treat lines in a list that includes 
"predicted" lines whose wavelengths are uncertain by ~10A (because one or both energy levels 
have not been measured in the laboratory), limiting the resolution of the flux distributions to 10 A. 
Moreover, they assume a solar iron abundance ~0.15 dex higher than the value currently accepted 
(Biemont et al. 1991; Holweger, Kock, & Bard 1995), and a solar ratio of light-element abundances 
[a/Fe] = 0. In contrast, metal-poor stars of the Galactic halo typically show [a/Fe] = 0.3 (Wheeler, 
Sneden, & Truran 1989), a factor-of-two enhancement in Mg, Si, Ca, and Ti. This affects the fit in 
the mid-UV, since it includes the very strong resonance doublet of Mg ii at 2800 A and a bound-free 
opacity edge of Mgi at 2512 A, plus many weaker Mg, Si, and Ti lines. 



2. Overview of This Work 

In this work we calculate the mid-UV spectrum line-by-line at high resolution, independent of 
the Kurucz flux distributions. Our line list does not include predicted lines, but rather is based 
on an updated laboratory line list augmented by manual addition of missing lines where necessary. 
To ensure the reliability of these calculations, we match as best we can the individual line blends 
observed in the mid-UV echelle spectra of real near-turnoff stars, spanning a wide metallicity range 
and including both solar and enhanced light-element ratios. Our purpose is to reproduce mid-UV 
spectra of turnoff stars accurately enough to determine T^s and to predict the composite mid- 
UV spectra of old stellar systems over their entire range of overall metallicity and light-element 
abundance ratio. 

In this first step, we describe this calibration procedure. Because the calibrations depend 
significantly on the stellar temperatures adopted, we have carried out new T^fi determinations in 
seven metal-poor turnoff stars and one mildly metal-poor hotter star. They are based on simulta- 
neously matching both mid-UV spectra (2280 - 3120 A) and the Ha profile and the strengths of 
high-excitation atomic lines in optical spectra. We thus do not rely on any photometry or parallax 
information in the T^g determination. Iron abundances and [a/Fe] ratios were adopted as derived 
from both weak and strong lines in the optical spectra. To assist in filling in the mid-UV line list, 
we have included mid-UV analyses of two stars with mid-UV spectra of very high resolution, the 
Sun and HD 128620 (= a Gen A), a well-studied nearby southern star with a greater-than-solar 
metallicity (Neuforge-Verheecke &: Magain 1997), for which T^fi was taken from the literature. 

In fitting each mid-UV spectrum, we adopted a single continuum normalization constant, the 
scale factor by which the observed mid-UV spectral flux was multiplied to match the flux predicted 
by theoretical calculation. This constant proved to agree to an average of 6% ± 5% with that 
expected from the ratio of stellar distance (found from the parallax) to stellar radius (found from 
the model log g assuming a reasonable stellar mass). Model U — B , B — V , and b — y colors all 
agreed extremely well with those observed, confirming Tpff errors of <50K. 

Our mid-UV and optical echelle spectra are described in §3. As discussed in §4, T^s was initially 
established from the optical spectra, and compared to the mid-UV fit. For the most metal-poor 
stars, this fit did not require the addition of missing lines, but did require the use of Kurucz models 
in which convective overshoot was turned off, in contrast to the Kurucz models in general use. 
To reproduce mid-UV line strengths and to match the pseudocontinuum of the mid-UV spectra 
of more metal-rich stars, mid-UV atomic line parameters were changed and missing lines were 
added. Transitions especially sensitive to the surface temperature were still poorly reproduced; this 
mismatch was reduced by adopting a mild enhancement of temperatures in shallow photospheric 
layers of the models of one-tenth solar metallicity and higher. 

The degree to which the resulting calculations match both mid-UV and optical spectra of 
turnoff stars is shown and discussed in §5. The match is very good for 2900 - 3120 A at all 
metallicities, and generally good throughout the mid-UV for metallicities at and below one-third 



solar. It could be improved at solar metallicities by obtaining additional high-resolution mid-UV 
spectra of a few well-chosen standards. 

Our Tcff results are presented in §6, where comparisons are made of model fluxes and colors with 
observations. We conclude in §7 with a summary and a discussion of the prospects for constraining 
age and metallicity of quiescent stellar systems by comparing such calculations to spectra of their 
integrated light. 



3. Observations 

We begin by listing the properties of the stars considered in Table 1. The first five columns 
following the HD number give the observed V magnitude, colors, and the parallax from SIMBAD. 
The next five columns, discussed in §6, show comparisons of observed properties of each star with 
those calculated from the models. The parameters of the model used to calculate both mid-UV 
and optical spectra of each star are given in the next four columns. In all cases but HD 128620 = 
a Cen A and the Sun, all the model parameters were redetermined here as described below. The 
final columns show T^s values determined previously in several works. 

Of the eight stars whose temperatures have been redetermined here from optical and UV 
spectra, seven are metal-poor main-sequence turnoff (MSTO) stars. The eighth, HD 128167, is 
a mildly-metal poor star lying above the halo main-sequence turnoff, either younger or a blue 
straggler. It was included to give temperature leverage to the mid-UV comparisons, as discussed 
below. Likewise the Sun and a Cen A were included to represent the high-metallicity end. The 
solar parameters were considered so well determined that no re-evaluation was necessary. The star 
a Cen A was not reanalyzed because no optical spectrum was available: it is inaccessible from the 
northern hemisphere. Its Tgfj and abundance have been determined to moderate accuracy from 
extensive abundance analyses, and its gravity as log g = 4.33 ± 0.02 as found from its visual and 
radial-velocity motion about its companion (Chmielewski et al. 1992; Neuforge-Verheecke &: Magain 
1997; Pourbaix et al. 1999). The iron abundance listed in Table 1 was adopted to best match its 
mid-UV spectrum. 

Optical spectra were obtained for all stars in Table 1 but the Sun and a Cen A . The Lick 
Observatory Shane 3m echelle spectrograph (Vogt 1987) was used with a TEK 2048 x 2048 CCD, 
at a FWHM resolution of 38,000. The spectra of HD 84937, HD 94028, and HD 106516 were kindly 
provided by Fulbright (2000). Exposure times on the 3m were typically lOmin for all stars except 
HD 128167 (40 s). Complete spectral coverage over 3900 A - 7800 A and S/N > lOO/pixel were 
obtained in all cases but one, although fringing compromises the accuracy of the data at the long- 
wavelength end of the range. For the one exception, HD 114762, the blue spectrum was obtained 
in two one-hour exposures by Tony Misch at Lick with the coude auxiliary telescope, using the 
same echelle spectrograph configuration. For the Sun, the flux spectrum of Kurucz et al. (1984) 
was downloaded from the Kurucz web site at http://cfaku5.harvard.edu; the full-resolution version 



was matched at 150,000 resolution. 

The mid-UV spectral observations of the Sun modeled here are the rocket spectra of Kurucz 
& Avrett (1981), recorded on film. Both passes of the center of the solar disk were used. Gaussian 
broadenings of 1.5 kms~^ for macroturbulence and 80,000 FWHM for other broadening sources 
provided a good match. A constant normalization scale factor of 1.10 was adopted, in keeping with 
the stated 10% uncertainty of its flux calibration. 

The stellar UV spectra were obtained with the Space Telescope Imaging Spectrograph (STIS) 
of Hubble Space Telescope (HST) from mid-1998 to mid-1999. Salient features are noted here; full 
details may be found from the HST archive at http://archive.stsci.edu. 

To obtain the mid-UV spectra of most stars, the near-UV MAMA detector was used with the 
STIS E230M grating over its default range of 2280 A - 3120 A at a FWHM resolution ~25,000. 
HD 106516 (= HR4657) was observed by Heap et al. (1998) in program #7433 with the 0.2 x 0.06 
aperture. Mid-UV spectra for HD 19445, HD 84937, HD 94028, HD 114762, HD 184499, and 
HD 201891 were acquired in snapshot program #7402 by Peterson &: Schrijver (2000), with two 
exposures obtained for HD 19445, HD 84937, and HD 114762. The 0.2 x 0.2 aperture was used and 
exposure times were 6 - lOmin. Consequently, S/N is lower, but flux levels are less dependent on 
pointing. Fluxes good to 1% for these six stars are indicated by previous observations of Mg ii line 
profiles in each star, made using the HST Goddard High-Resolution Spectrograph (GHRS) and the 
2x2 aperture in program #5869 (Peterson & Schrijver 1997). 

Mid-UV spectra of the stars HD 128167 (=HR5447), from #7433, and HD 128620 (= HR5459 
=Q;CenA), from #7263 (Linsky et al. 2000), were recorded in multiple exposures with the E230H 
grating. Spectral coverage is 2380 A - 3160 A for HD 128167 and <2280 A - 3160 A for HD 128620. 
A FWHM resolution of ~60,000 was achieved with the 0.2 x 0.06 aperture for HD 128167, and the 
0.2 X 0.05ND aperture for HD 128620. For neither star do fluxes always reproduce exactly in orders 
obtained in more than one exposure; differences as large as 10% are seen. Similar uncertainties are 
anticipated for HD 106516, where the 0.2 x 0.06 aperture was also used. 

Rotation is detected in at least two stars and possibly in two more, from the extra broadening 
required to match line profiles. For HD 128167, HD 106516, HD 84937, and HD 114762, line breadths 
were matched in the spectral synthesis by adopting rotational velocities v sin i = 9, 7, 5, and 4 km s^"*^ 
respectively. Thus the resolution of the mid-UV spectrum of HD 128167 is effectively degraded to 
that of the stars observed at lower resolution, affecting the discernment of mid-UV lines. This star 
nonetheless has one of the sharpest-lined spectra among those obtained in program #7433. 



4. Spectrum Analysis 

Spectrum analysis was based on visually matching each observed spectrum to ab initio radiative- 
transfer calculations assuming static models in LTE. The program SYNTHE (Kurucz & Avrett 



1981), the Kurucz (1995) and Castelli et al. (1997) grids of ATLAS9 models, and the Kurucz hne 
hsts for atomic species and molecular hydrides were downloaded from the Kurucz web site, where 
further details on each of these can be found. As described below, all models, including that of the 
Sun, were ultimately drawn from the Castelli et al. (1997) grid alone. The more metal-rich ones 
were modified by introducing enhanced temperatures in the shallow layers. 

The SYNTHE program accepts as input both line lists and a model atmosphere with tempera- 
ture and other quantities tabulated for each depth, and calculates an array of spectra at individual 
emergent angles. These may then be used individually (as for spectra of the center of the solar 
disk) or convolved with a zero or positive rotational velocity to generate an emergent flux spectrum. 
In either case, the resulting spectrum can be broadened by Gaussians (or other profile shapes) to 
incorporate macroturbulence and instrumental broadening. 

SYNTHE was modified to run under UNIX on a Sun Ultra-30 at the University of Virginia. 
Mid-UV spectra from 2280 A to 3160 A were calculated at a resolution of 330,000 (500,000 for the 
Sun and HD 128620), and optical spectra at a resolution of 500,000. A Gaussian macroturbulent 
broadening of 1.5 kms~^ was assumed for all stars as well as for the Sun. Spectra for some stars 
were rotationally broadened as listed in the last paragraph of §3, as was the solar flux spectrum, 
for which a mean rotational velocity of 2 kms^^ was adopted. The spectra were then broadened 
by Gaussians corresponding to the values of instrumental resolution listed individually in §3. 

The line lists and the source code of the opacity routines available at the Kurucz web site give 
specifics on the wide variety of opacity sources treated by SYNTHE. For example, the significant 
absorption by the bound-free absorption edge of Mgi near 2512 A is treated both as line opacity, 
with the line list including the bound-bound transitions to A^ ~ 60, and as continuous opacity 
blueward of the bound-free edge at Uq = 39759.8 cm~^, by the subroutine MGIOP. Near Vq, this 
routine adopts a cross section of 40 Mbarn with a (z^/fo)~^^ frequency dependence, based on the 
laboratory measurement of 46 it 12 Mbarn by Lombardi, Smith, & Parkinson (1981) and the 
convergence of high-level lines. Where this formalism falls below the theoretical cross section, the 
program adopts the latter, with a cross section of 20 Mbarn and a {u/uo)~'^ frequency dependence 
(e.g., Butler et al. 1993). 

For calculations redward of 5000 A, the line list of Peterson, Dalle Ore, & Kurucz (1993) based 
on the solar spectrum was adopted, with modifications due to the use here of a flux, rather than 
central intensity, spectrum for the Sun and to the incorporation of an approximate chromosphere 
in the solar model. This led to the deduction of larger gf-values for the strongest lines except 
where laboratory gf-values were adopted, as for low-excitation Fe i lines with furnace measurements 
(Blackwell et al. 1982). 

Blueward of 5000 A, the Kurucz 25 May 1998 atomic line hsts a/0300.100, 5/0400.100, and 
5/00500.100 were used. These differ significantly from those used in the Kurucz flux distribution 
calculations in including only atomic lines identified in the laboratory (with laboratory gf-values 
where available) - like those of Kurucz & Bell (1995) but with the addition of many lines of Fe i, Fe ii. 



and other species whose upper energy levels were newly determined in near-UV laboratory spectra 
(e.g., Nave et al. 1994). Because the strengths or profiles of many lines were mismatched when 
using line parameters directly from this list, their gf- values and damping constants were changed 
where necessary, as described shortly. Since many lines appear in the STIS mid-UV spectra that 
were not modeled by lines in the list, such missing lines were added blueward of 3120 A. Redward 
of this, no missing lines were added, but gf- values were changed to match stronger lines in the Sun 
in the regions depicted in Fig. 4 below. 

The Tcfj, log g, [Fe/H], xt, and [a/Fe] values for each reanalyzed star were initially established 
from Ha and weak lines of neutral species in the optical spectra, reassessed from optical lines of 
once-ionized species, then confirmed or revised from the UV. For the two most metal-poor stars, 
mid-UV spectra provided an immediate check on Tcs . For the others, the following steps were 
iterated until all diagnostics converged. 

First the Ha profile and an assumed log g were used to estimate Tgfj. [Fe/H] was adjusted until 
weak Fei lines in the 5100 A - 5400 A region were well matched; their strengths are largely gravity- 
independent. The strength of weak Til, Mgi, and Sii lines determined [a/Fe] . Microturbulent 
velocity Xt was set by matching weak and strong iron and titanium lines. Gravity log g was then 
checked by demanding matches for Mgi, Mgii, and the wings of the Mg6 lines, for Fei and Fell, 
for Til and Till, and for Si I and Sill. This also constrains Tcfj , for the Mgll and Sill doublets 
arise from highly excited lower levels (>8eV), while the optical Tin and Fell lines originate from 
much lower levels (<3eV). Roughly two dozen Fei lines, a dozen Til lines, and a dozen total Tin 
and Fe II lines were detected for the most metal-poor stars, with double this number typical of the 
more metal-rich stars. The Mgii and Sill doublets at 4481 A and at 6347.1 and 6371.3 A were seen 
in all cases, but Sill only marginally so in HD 84937 and HD 19445. 

Once the diagnostics began to converge, iterative adjustments were made to the Kurucz mid- 
UV atomic-line parameters. First, gf-values were modified for atomic lines identified in the labora- 
tory. The relative contribution of identified lines to a blend was ascertained from the behavior of 
the blend in the hot star HD 128167 versus the cooler ones, and in the higher-gravity stars such as 
HD 19445 versus the lower-gravity ones. Damping constants were changed in some cases to match 
line profiles. Damping constants rather than gf-values were changed for the strong resonance lines 
well-observed in the interstellar medium for which gf-values are tabulated by Morton (1991). 

Next, for lines lacking identification, the laboratory line list was searched for transitions close 
enough in wavelength to perhaps be responsible. The gf-value of a nearby identified line was 
increased by up to 4.5 in the log, and the calculation repeated. This line was dropped if still not 
strong enough; otherwise its gf-value was adjusted to match. The high resolution of the a Cen A 
and solar spectra showed that most such identifications were very well matched in wavelength, and 
thus likely to be correct, when based on increases of < 2.5 in the log. 

For the two most metal-poor stars, these steps were sufficient to enable the continuum to 
be defined unambiguously throughout most of the mid-UV region. Its definition in more metal- 
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rich stars was hampered by hne crowding, yielding a pseudo continuum that the modehng did not 
properly reproduce. 

Much of this pseudo continuum mismatch appeared to be due to individual transitions missing 
from the line list. This was indicated by the very common occurrence throughout the mid-UV of 
unmodeled absorption features found at the same wavelengths in all stars, features whose strength 
steadily increased in going from weak-lined to strong-lined spectra. Most of these transitions are 
weak, consistent with the fact that the intrinsically strongest transitions of each species are the 
ones most likely to have been detected in the laboratory. 

Except just redward of the Mgi edge at 2512 A, most missing lines are probably due to iron. 
Its abundance is high and the spectra are rich in lines of Fe i and Fe ii, which are well-populated 
species in solar-temperature stars. Of the elements with similar first ionization potentials, only 
magnesium and silicon have comparable abundances; the mid-UV spectra of their low-ionization 
stages are generally much less complex, as are those of the more abundant CNO elements. While 
iron-peak elements such as Cr, Mn, and Ni have many transitions throughout the mid-UV, iron 
dominates because it is twenty or more times as abundant. In most halo stars above one-tenth 
solar metallicity, the relative abundances of these elements with respect to iron are within 0.2 
dex of those of the Sun (although larger differences may occur, especially at very low metallicity: 
McWilham et al. 1995; Nissen et al. 2000; Prochaska k McWilham 2000). The mismatch that wih 
result from assigning missing lines of these elements to iron should be mild. 

In the region immediately redward of the bound-free Mgi opacity edge at 2512 A, the proba- 
bility is significant that missing lines may be due to magnesium rather than iron. Redward of the 
edge there are a very large number of transitions to known highly excited levels of the Mg i atom, 
and presumably a very large number of similar transitions to closely-spaced highly excited levels 
that are as yet unidentified in the laboratory. If many weak transitions redward of 2512 A are 
due to magnesium but have been assigned to iron instead, our calculations will underestimate the 
absorption in this region in metal-rich alpha-enhanced populations. We see below that this region 
is poorly modeled in any case once metallicity exceeds one-tenth solar. 

The modeling then proceeded iteratively by assigning tentative identifications and transition 
probabilities to each missing feature, then recalculating until an acceptable match was achieved or 
deemed unattainable with existing data (as for the regions just redward of 2550 A and 2640 A in the 
solar spectrum). Features were attributed to either low-excitation Fei, moderate-excitation Fei, or 
Fe II, whose line strengths depend differently on stellar temperature and gravity (pressure) because 
of the Boltzmann and Saha equations. With a Tcs ~1000K higher than those of other moderately 
weak-lined stars, HD 128167 proved exceptionally useful in making these assignments. The higher 
resolution of the spectra of the strong-lined stars HD 128620 = a Cen A and the Sun was essential 
in discerning the wavelengths of the many features emerging above one-third solar metallicity. 

A brief examination was made of the gf-value changes for low-excitation Fe i lines, namely 
those identified in the laboratory with (air) wavelengths falling between 2600 A and 3000 A, and 
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with a lower level of 6000 - 10000 cm^^ and an upper level of 40000 - 50000 cm^^. Sixty-five such 
Fei lines have gf- values from the laboratory measurements of O'Brian et al. (1991); 52 were not 
changed, while 13 were changed by from —0.6 to +0.4 dex. Nine lines have gf-values from the 
laboratory measurements of Fuhr, Martin, & Wiese (1988); five were unchanged, and four changed 
by from —0.7 to —0.3 dex. Thus the mid-UV gf-values found in this work generally concur with 
the mid-UV laboratory measurements. 

In contrast, fifteen lines have theoretical gf-values from Kurucz (1994); four were not changed, 
but 11 were changed by from —2.2 to +0.6 dex, all but one decreased. The theoretical gf-values 
of Kurucz (1994) are seen to be subject occasionally to very large errors, usually overestimates. 
This can be attributed to the finite size of the Cray computer used in those calculations, which 
dictated that the number of levels included be truncated, and so upper levels were preferentially 
omitted from the model atom. Because the sum of the gf-values of transitions arising from a 
given lower state is normalized to unity, gf-values for the transitions included in the calculation 
will be overestimated because of the exclusion of others arising from the same level. Thus the 
calculations will overestimate the gf-value by increasing amounts as the fraction decreases of the 
actual possible upper states that are included in the atomic model. This fraction decreases to 
the blue, because higher upper levels are involved for a given lower level. Indeed, the severity of 
overestimates appears to increase to the blue: all the above changes made to the Kurucz (1994) 
gf-values were —1.5 or lower below 2650 A, but redward of this only one change this negative was 
required. The overestimates among the Kurucz (1994) Fei gf-value calculations are thus seen to be 
due to incompleteness and not to more fundamental difficulties. Undoubtedly many of the weaker 
mid-UV transitions are missing entirely from the predicted list. 

For the hydrides, separate treatment was required, as the lines are naturally associated into 
bands, and the Kurucz hydride calculations are now several decades old. For OH, the theoretical 
calculations of Gillis et al. (2001) (communicated in August 2000 by M. Bessell) were used wherever 
possible. For CH, the LIFBASE data (e.g., Luque & Crosley 1999) were used (communicated 
in August 2000 by M. Bessell with lines cross-identified). Since these cover only the strongest 
bands, calculations must still rely on the Kurucz values for the others, notably the 4-3 band of 
OH with many weak but discernible lines in the Sun. A comparison of the Kurucz OH transition 
probabilities against those of Gillis et al. (2001) showed strong trends, with discrepancies increasing 
at high rotation number J. To place them on approximately the same scale, the Kurucz transition 
probabilities for the OH and CH lines were decreased by an amount 6 = —0.25 — 0.05 * Ji dex. Once 
missing atomic lines were largely filled in, a comparison was feasible of calculated OH and CH 
line strengths against those of the Sun and stars. Because solar OH lines were otherwise modeled 
too strong, all the transition probabilities were reduced by an additional 0.15 dex uniformly. In 
part because weak OH lines in the Sun were reasonably well modeled but strong ones were still 
too strong, a chromospheric mitigation of the surface temperature drop was invoked as described 
shortly. As illustrated below in Figs. 3d and e, these changes suffice to reproduce both strong and 
weak 0-0 OH features in stars as well as the Sun. 
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The calculations were first run using Kurucz (1995) models, but the very metal- poor stars 
showed a Tefj deduced from Ha as much as 200 K higher than that from fitting the mid-UV. That 
discrepancy was removed by adopting Castelli et al. (1997) models, models in which convective 
overshoot has been turned off. We did not use the alpha-enhanced Castelli et al. (1997) models, as 
they are available for only a single value of the alpha-element abundance enhancement. 

Like the Kurucz (1995) models, the Castelli et al. (1997) models we used were calculated with 
Xt = 2 kms~^ and a logarithmic solar iron abundance of —4.37 with respect to the abundance 
by number of hydrogen plus helium, which amounts to —4.33 with respect to hydrogen alone. 
These choices are about 1 kms~^ higher in xt than the ~1 kms^^ typically found for the solar 
flux spectrum, and about 0.15 dex higher in iron abundance than indicated by the meteoric iron 
abundance and by several recent determinations for the Sun (e.g., Biemont et al. 1991; Holweger, 
Kock, & Bard 1995). Thus both sets of models were calculated with higher line blanketing than 
expected for a solar-type star of the nominal model metallicity, since iron lines dominate the UV 
and optical spectrum, and the higher xt increases the strengths of moderately strong lines of every 
species. 

For each star including the Sun, a Castelli et al. (1997) model was interpolated to the nearest 
50 K in Tcff (25 K for the Sun), 0.1 dex in log g^ and 0.05 dex in [Fe/H]. The microturbulent velocity 
Xt was changed at all depths to the value indicated for each star in the table and figures below. To 
minimize the overblanketing effect, our models were interpolated to a nominal grid abundance 0.10 
- 0.15 dex lower than that desired. To set abundances, the abundance of the interpolated model 
was then raised to the desired value and the relative abundance of iron dropped by 0.15 dex. The 
specific light-element abundance elevations required by the optical observation were adopted, or 
assumed to be zero for the Sun and a Cen A , and all species equilibria were recalculated. This 
procedure incorporates the specific Xt and iron and light-element abundances set for the model 
in calculating both the ionization and molecular equilibria and the line and continuous opacities 
for its spectrum. Concerning the model structure, the procedure compensates for the effects of 
overblanketing due to the high choice of iron abundance, but not for the light-element ratio nor 
for the microturbulence. We anticipate very mild errors in model colors and temperature structure 
until spectra become dominated by molecular lines or otherwise are considerably stronger-lined 
than any of these. The comparison of model and observed colors (§6) bears out this expectation. 

For the Sun, a model with T^s = 5775 K, log g = AA, [Fe/H] = 0, [a/Fe] = 0, and xt = 
1.0 kms^^, was interpolated in the above way. For HD 128620 = aCenA we adopted a similarly 
interpolated model with Teg = 5800 K, log g = 4.3, [Fe/H] = -^0.15 dex, [a/Fe] = 0, Xt = 1-0 kms"^ 
. While Neuforge-Verheecke & Magain (1997) found T^s = 5830 K and [Fe/H] = 0.25 dex, most 
other analyses have derived somewhat lower T^g and [Fe/H] values. At T^s = 5800 K and xt = 1-0 
kms~^, [Fe/H] = 0.15 dex provides a better match to the mid-UV spectrum of a Cen A. 

Even with the use of Castelli et al. (1997) models interpolated as above, another major type 
of mismatch persisted in all strong-lined stars: cores of lines formed near the surface appeared 
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too strong in the calculations, especially in extremely strong mid-UV lines such as those of Mgll 
and Mgi or Fell and Fei. This is definitely not due to wayward gf-values or damping constants, 
because of the trend seen in the mismatch. Profiles calculated for the strongest lines in the most 
metal-poor spectra generally matched or were somewhat too high in the core, but were somewhat 
too low in the core in stars of intermediate line strength, and became too low across the the central 
ten Angstroms or more in the strongest-lined stars. (See Fig. 3a below for residual trends of this 
nature.) As noted above, a related mismatch was seen within the OH spectrum itself: in any 
given strong-lined spectrum, notably that of the Sun for which the model parameters are extremely 
well known, weak OH lines were matched at a higher oxygen abundance than strong lines. The 
possibility discussed in §5 that the solar continuum may be drawn too low would lead to a mismatch 
of opposite sign. To the extent that the gf-values are reliable, then, this indicates an error in the 
temperature structure of the model. 

The most obvious possibility is a chromosphere, which is not incorporated in flux-constant, 
static models such as those of SYNTHE. The solar chromosphere, for example, is seen in images of 
the limb and evident from emission reversals in the cores of strong spectral lines such as Mg II. It is 
generally modeled by an increase in temperature in relatively shallow layers of the atmosphere (e.g., 
Fontenla et al. 1999). Any of several non-radiative heating sources extending from the photosphere 
through the chromosphere into the corona could be responsible (e.g., Schrijver 1995; Schrijver et al. 
1999; Sterling 2000). 

We consequently changed the temperature of individual points in each interpolated stellar 
model with [Fe/H] > —1.0, beginning just below the point where the model becomes optically thin 
at 5000 A (typically at depth number 51 or 52). (We tried this for the more metal-poor models as 
well, but the profiles of all the strong lines became much too high near the core, so we left these 
models as they were.) To best fit both optical and mid-UV spectra simultaneously, we increased 
the temperature at model depths 28 through 48 by 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 110, 
120, 130, 140, 150, 150, 140, 120, 90, 60, and 30 K respectively. This temperature enhancement 
may be considered an approximate chromosphere, though it is based entirely on empirical fits of 
one-dimensional spectral line calculations to observed spectral line profiles. It incorporates neither 
energy transport considerations nor matches to spatially resolved structures, both of which lie 
beyond the scope of this work. 



5. Matching Mid-UV Spectral Calculations 

The degree to which our mid-UV and optical spectral calculations match spectra of real stars 
is illustrated in the figures. The first three show mid-UV spectra, and the last, optical. Within a 
figure, panels are ordered by wavelength. Flux (central intensity for the solar mid-UV spectrum) 
is plotted versus wavelength in air. The spectrum calculated for each star (light line) is plotted on 
top of that observed (heavy line) . 
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In Figure 1, the global mid-UV fit is shown for the nine stars with absolute fluxes. All spectra 
have been Gaussian smoothed to a FWHM resolution of 500 (550 for HD 128127 and aCenA). 
Each spectrum is displaced vertically by -^20% from the one below. For each star a single mid-UV 
continuum normalization constant was adopted, as listed on the left. On the right, the HD number 
of each star appears in bold. Below are found the model-atmosphere parameters adopted for the 
corresponding calculation: T^s (in K), C wherever the chromosphere was modified as described at 
the end of §4, log g, [Fe/H], xt (in kms~^). On the following line appear the values in dex of the 
light-element abundance enhancements [Element /Fe] (in dex), first that of oxygen, then those for 
Mg, Si, Ca, and Ti. (The latter four are grouped together wherever they are the same.) While the 
fits are generally good, serious discrepancies are noted in a Cen A at the bottom of the figure near 
2550 - 2600 A and 2640 - 2710 A, regions of relatively little line absorption in the more metal-poor 
spectra. 

Figure 2 illustrates the effect on mid-UV fits when model parameters are changed somewhat. 
Two of the stars of Fig. 1 are depicted: the moderately weak-lined HD 94028 (top two spectra) and 
the rather strong-lined HD 184499 (bottom three). Star and model identifications and normalization 
constants are indicated as in Figure 1. The lower plot for each star represents its best fit, a replot 
of that star's comparison in Figure 1. 

Comparing the bottom plot to the middle plot shows the effect for HD 184499 of decreasing 
[Mg/Fe] by 0.3 dex to the solar [Mg/Fe] ratio. Although minimized by the adoption of a normal- 
ization constant 11% higher, the mismatch is still evident. The strengths of the wings of the strong 
Mgii doublet at 2800 A (and of Mgi at 2852 A) are poorly reproduced, as is the overall pseudo- 
continuum level near and below the bound-free edge of Mgi at 2512 A. Similar calculations run 
for HD 94028 show a more dramatic effect in the Mgii and Mgi line wings, which are weak enough 
that little renormalization is possible, but a smaller effect near 2512 A, where the H^ opacity is 
dominant because of the reduced magnesium abundance. 

Comparing the plot just above the best fit for each star to the best fit itself shows the degeneracy 
of the mid-UV spectra with respect to a simultaneous change in Tcs , log g and [Fe/H] designed 
to preserve overall line strength and ionization equilibrium. For each plot just above the best fit, 
Tgff was decreased by 100 K, [Fe/H] and log g both concurrently decreased by 0.1 dex, and the 
spectrum renormalized accordingly. This change is very difficult to detect from fitting the mid-UV 
slope alone at these metallicities and this resolution. The increased metallicity reduces the blue 
pseudocontinuum relative to the red at the same time that the increased T^g raises it. The change 
is marginally detectable in HD 94028 from the mismatch of low-opacity windows, but not in HD 
184499 because the overall match is less satisfactory. Such a change does lead to a reduction in the 
normalization parameter, by 11% for HD 94028 and 17% for HD 184499. In principle this might 
be detected by other means, but only if [Mg/Fe] is properly modeled. For as seen above, taking 
[Mg/Fe] = (as in the Kurucz flux distributions) can largely compensate in stars such as HD 
184499 that show a light-element enhancement. 
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Our conclusion, then, is that if the niid-UV norniahzation constant is treated as a completely 
free parameter (as here), the mid-UV flux distribution by itself can currently be used to set T^fi only 
in MSTO stars with [Fe/H] < —1.5. While the mid-UV was useful in demonstrating the need for 
models with no convective overshoot and with an approximate chromosphere, our determinations 
of Tefj for the more metal-rich MSTO stars have relied on the optical spectra instead. We show in 
Fig. 4 below that the Balmer line profiles and the strengths of other high-excitation lines are well 
matched in the optical spectra of all stars at the same Tefj as the best-fit mid-UV spectrum, and in 
§5 that model colors match observations as well, so that T^s is indeed very well established here. 
However, our goal of T^q determination for turnoff stars of all metallicities from mid-UV spectra 
alone is not yet met. We are not able to calculate the low-opacity windows blueward of 2900 A 
accurately enough at near-solar metallicity, as we now demonstrate. 

Figure 3 plots at higher resolution and higher scale the same mid-UV spectral comparisons as 
are shown in Figure 1, plus that of the Sun. Panel a shows 2627.5 - 2630.5 A, dominated by strong 
lines. The other panels show regions of locally low line blanketing: b, 2642.5 - 2645.5 A; c, 2900.0 
- 2903.0 A; d, 3082.5 - 3085.5 A; e, 3100.0 - 3103.0 A. In all plots in Fig. 3, the spectrum of HD 
184499 is offset by 30% from that of HD 128620 = a Gen A, which in turn is offset by 50% from 
that of the Sun, for which 1.0 marks the height of the true continuum. The same normalization 
constant as in Fig. 1 has been adopted for each stellar spectrum, as was a normalization for the 
Sun 1.10 times that given in the Kurucz & Avrett (1981) atlas, as noted on the left. Identifications 
are provided at the top for the strongest lines as calculated in the solar spectrum. Manually added 
lines are flagged with a colon following the decimal digits of the wavelength in Angstroms. Next 
appears the species, the lower excitation potential in eV of atomic transitions or the band and line 
of molecular ones, the digits following the decimal of the residual flux in the unbroadened line of 
the solar spectrum calculation, and the log gf-value adopted. 

In Fig. 3e of the 3100.0 ~ 3103.0 A region, calculated and observed spectra agree reasonably 
well in all stars, including aCenA and the Sun. No "missing" lines are strong enough to merit 
explicit identification at the top of the figure. True continuum is reached in each spectrum, but 
fleetingly, near 3102.5 A; its placement is sensitive to the treatment of adjacent weak lines. An 
unblended 0-0 OH line appears at 3101. 230 A, next to a near- continuum high point; both are 
modeled well in all spectra. The strong atomic lines are well reproduced, as are two moderately 
strong CH lines at 3100.120 and 3100.170 A. 

Fig. 3d offers a further picture of how well both continuum and OH features are modeled. 
Another continuum window appears near 3084.7 A, but again with a weak line superimposed in 
the solar and aCenA spectra. The relatively unblended OH 0-0 band line at 3084.896 A, used 
in several oxygen-abundance investigations, is well reproduced over the factor-of-100 metallicity 
range of the stars in whose spectra it is detected, including the Sun itself. Nearby at 3085.2 A is a 
blend dominated by a stronger 0-0 OH line; it too is reasonably well matched in all stars. Recall 
that to achieve this agreement, the Gillis et al. (2001) theoretical gf-values for OH were reduced 
by 0.15 dex; otherwise these features and that of Fig. 2e, plus the overwhelming majority of other 
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OH lines, would be too strong in the Sun. That is, unless the solar continuum were raised. 

The solar continuum appears to be well defined, and the OH lines of all bands matched to 
~0.1dex, throughout the 3000 - 3100 A region. Because there are so few continuum windows, 
however, we cannot rule out a systematic underestimate of the continuum in a Cen A and the Sun 
of perhaps 5% or 10%. Roughly speaking, raising the a Cen A mid-UV continuum would allow it to 
be modeled at an [Fe/H] value more consistent with that of Neuforge-Verheecke &: Magain (1997), 
while raising the solar continuum might eliminate the present need to reduce the theoretical OH 
gf-values by 0.15 dex. As noted at the end of §6, it would also bring the mid-UV normalization 
constants of the Sun and a Cen A closer to their expected values. A definite answer could be 
obtained by deriving the stellar parameters of a Cen A from an optical spectrum in the same way 
as for the other stars, and from modeling an improved spectrum of the solar near-UV. 

Random continuum variations approaching 5% occur for two, and possibly all three, of the 
stars that were noted in §3 as observed through a smaller aperture. For the continuum seems 
appropriately placed for all stars in Fig. 3e, but in Fig. 3d the continuum of the HD 106516 observed 
spectrum appears too low, and in Fig. 3c that of HD 128167 may also be too low. Such variations 
might also be present in the HD 128620 = a Cen A spectrum as well, but remain undetected because 
of the difficulty of matching its multitude of lines. 

Fig. 3c shows another low-opacity window at 2900 A . Here the number of missing lines added 
remains small, and agreement is nearly as good. OH lines are less important, being generally minor 
constituents of blends. However, windows of true continuum have disappeared from the solar and 
a Cen A spectra, and are fleeting at best in the echelle spectra of lower resolution that were acquired 
for mildly metal-poor stars. Continuum placement depends not so much on OH features but rather 
on whether all missing lines have been spotted. Clearly that is not the case at solar metallicity, 
where weak missing lines apparently overlap to the extent that they cannot be distinguished as 
individuals. Although the cumulative effect of their neglect on the total absorption is small, is it 
significant for the continuum placement, which consequently cannot be determined as well as was 
the case above from fits such as this. It is placed here and for the remaining panels by adopting 
the same normalization factors required to match Figs. Zd and e. 

By 2645 A, in Fig. 36, the fits have deteriorated considerably. This is especially true at and 
above solar metallicity, where the numbers of missing and identified lines noted at the top are 
comparable. True continuum still appears in the weakest-lined spectra near the top of the panel, 
but elsewhere it is gone, lost to a forest of weak lines. Although the fit to the metal-poor stellar 
spectra is still reasonable, the high-resolution spectra of the a Cen A and the Sun reveal that too 
few missing lines have been added, often too strong and in the wrong places. For the adjacent 
region from 2645 to 2700 A, the sheer number of missing lines precluded their reliable identification 
at the low resolution of the low-metallicity stellar spectra. Because many were simply omitted, the 
calculations overestimate the spectral flux in that region - mildly in moderately metal-poor stars 
but seriously at and above solar metallicity, according to Fig. 1. 
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From Fig. 3a, it is clear that these difficulties are reduced in regions dominated by strong 
atomic lines. There the fit is dictated more by the ability of the calculations to match the cores and 
the damping wings of strong lines, and not by modeling the weak or missing features. These are 
largely suppressed: although the region shown is blueward of that of panel b, only one missing line 
is noted in it. The wings of the strong lines are reproduced well in the four top spectra but begin 
to be somewhat too broad near the core in HD 106516 and HD 201891. The entire itlA about 
the core is too deep in a Cen A and the Sun. This could perhaps be remedied by continuing the 
temperature enhancement to shallower surface layers than those indicated in §4. 

Figure 4 shows fits to the optical spectra of the Sun and all stars except a Cen A. Identifications 
are given for the strongest lines in the solar spectrum. Since the optical spectra were not fluxed, the 
continuum was normalized to unity. The same models were used to calculate the synthetic spectra 
as in Figures 1 and 2. The wavelength regions in panels a - e depict a variety of magnesium lines 
and those of other heavier elements, to illustrate the extent to which the spectral calculations match 
features of a wide range of excitation, ionization, and strength. 

Figure 4/ shows Ha , to illustrate the satisfactory choice of T^s . The normalization of the Ha 
profile was aided by its position near the center of an echelle order, and was done by interpolating the 
normalization functions of the two adjacent orders, centered ±75 A away. This was straightforward, 
for the Ha line is much narrower than this and the continuum is very well defined by a spline fit 
of low order; moreover, the normalization functions for the adjacent orders were very similar. Note 
that most absorption lines in this panel are telluric, not stellar, as evidenced by their shift in 
apparent wavelength from one star to the next. 

Included in Fig. 4/ above the best fit for HD 94028 is the spectrum for the cooler HD 94028 
model whose mid-UV spectrum is shown in Fig. 2. This illustrates how the fit to Ha deteriorates 
when models of 100 K lower Tcfj are adopted. The Ha profile is relatively insensitive to gravity 
in this case: calculations adopting a model that differed from the best-fit model only in having 
log 5f 0.2 dex larger produced a profile that deviated only near 6560 A and 6565 A, by about 
the thickness of the heavy line. Similar calculations showed no sensitivity of the Ha profile to the 
chromospheric elevation of temperature, nor to changes of 0.3 dex in iron abundance or light-element 
ratio. However, Kurucz (1995) models required a Tefj 200 K hotter, thanks to their inclusion of 
convective overshoot. 

Panels a - e of Fig. 4 show that all the optical magnesium features are well matched using the 
same magnesium abundance and the same model as for the best fit to the mid-UV and to Ha. Panel 
a plots the region near the very high-excitation Mgii doublet at 4481.2 A; b, the ground-state Mgi 
line at 4571.1 A; c and d, the Mg6 lines at 5167.3, 5172.7, and 5183.6 A; and e, the high-excitation 
Mgi line at 5711.1 A. These panels also demonstrate that the calculations match strengths and 
profiles of other atomic lines regardless of species, excitation, or strength. Among the lines depicted 
are Til at 5173.7 A, Tin at 5185.9 A, ground-state Fel at 5166.3 A, Fell at 4576.3 A, and Sii at 
5708.4 A. These panels also contain missing features: not modeled at all in the calculations are the 
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absorption lines seen near 4480.8 A, 5170.7 A, and 5181.3 A. They are not telluric, since they appear 
at the same position in the stellar spectra, and are especially visible in the two stronger-lined stars. 



6. Results and Comparisons 

Our results for the parameters of the eight metal-poor turnoff stars analyzed here are listed 
in columns 12 - 15 of Table 1. The T^g values found for these stars by works cited in §1 are 
included for reference. Our Tefj values are seen to fall in the middle of the range of previous results, 
preferring neither the hotter nor the cooler values. Unfortunately, the number of stars is too small 
and the scatter too large for more definitive conclusions. 

The reliability of our T^s values is supported by a comparison of model colors to those observed. 
For each star including the Sun, we interpolated the UBV and h — y colors and V magnitudes 
calculated by Castelli et al. (1997) to the Tcfj, log g, and [Fe/H] values of each model we used. We 
have made no reddening corrections, nor any corrections for the temperature enhancement invoked 
in the more metal-rich models as described at the end of §4. The ninth, tenth, and eleventh columns 
of Table 1 show the difference m. B — V , U — B , and b — y between the observed and calculated 
colors for each star. The agreement is extremely good in all colors over the entire metallicity range. 
The mean difference is —0.010, —0.018, and 0.000 mag in B — V , U — B , and b — y respectively. 
The standard deviation of an individual difference is 0.010, 0.026, and 0.015 mag, marginally larger 
than typical observational errors alone. The color differences are small even for aCenA, for which 
both metallicity and Tefj are influential (we interpolated to a model [Fe/H] = -1-0.05), and for which 
we have no optical spectrum and so did not derive parameters. For the other stars, B — V and 
b — y are sensitive primarily to Tefj; their deviations correspond to ±45 K with no allowance for 
observational error. Their observed colors are slightly bluer than those of models, but the size of 
the differences points to a scale error in Tefj of less than 50 K . 

We thus estimate an uncertainty of ±50 K in each Tgff determination, from this and the poorer 
quality evident in Figs. 2 and 4/ of the fits to mid-UV fluxes and Ha profiles as Tefj is increased by 
100 K, plus the deterioration of the fit shown in Fig. 4a to Mgii. From plots like those of Fig. 4 at 
other wavelengths, we estimate the uncertainty in [Fe/H] to be ±0.05 dex at fixed Tefj , increasing 
to ±0.1 dex when the Tcs uncertainty is included. The mid-UV slopes and Ha profiles are only 
weakly sensitive to log g, so our results for log g are susceptible (as are those of other groups relying 
on ionization equilibria) to errors in gf-value scales and to possible non-LTE effects. We have not 
yet taken full advantage of the technique of matching damping wings of strong lines, incorporating 
improved treatments of hydrogen broadening such as those of Anstee &: O'Mara (1995) and Barklem 
&; O'Mara (1997). It would be valuable to do so given the constraints this would offer on alternative 
measurements of the fundamental properties of these standard stars, for example parallaxes from 
HIPPARCOS and other space missions. However, the following comparisons of observed and model 
fluxes indicate that the log g values are good to ±0.1 dex as they stand. 
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We first compare the UV normalization constant with that expected. According to Mihalas 
& Binney (1981), the ratio of observed to emitted stehar flux is f /F = 0^/4, where the angular 
diameter 9 = 2R/D^ R being the stellar radius and D its distance. For each star we calculate 
log R in solar units as 0.5(log M — log g + 4.4377), assuming masses M = 0.8 Mq for for the metal- 
poor stars based on stellar isochrones (VandenBerg et al. 2000) and 1.15 Mq for HD 128167, and 
adopting 1.16 Mq for aCenA as derived from its orbit by Pourbaix et al. (1999). The distance in 
solar units is found from the parallax, and the theoretical flux normalization constant follows from 
F/ f times 4.2545 xlO^*^ to convert to square steradians. The ratio of this theoretical constant to 
that actually used to scale the observations in Figs. 1 and 3 is listed in column 7 of Table 1. 

Agreement is good for all stars but one. For aCenA, the ratio is 0.46, suggesting an observed 
flux that is too low. Indeed, Linsky & Wood (1996) saw this in their 1995 mid-UV aCenA HST 
observations with GHRS: relative to earlier GHRS and lUE observations of a Cen A, "our fluxes are 
lower (by factors of 5.6 and 4.8 for the Mgii and Fell spectra, respectively), presumably because 
of the star's being misplaced at the edge of the small aperture." From a visual comparison with 
Fig. 1 of Linsky &: Wood (1996), the fluxes we are using for the Linsky et al. (2000) STIS a Cen A 
spectrum should be raised by a factor of 2.1. This brings the a Cen A ratio to 0.96, and results in 
an average of the nine stellar ratios of 1.06 it 0.05, with an individual standard deviation of 0.15. 
For the Sun, the ratio is 0.91 given the scale factor of 1.10 noted above, which is within the stated 
normalization uncertainty but still low. These a Cen A and solar ratios would be closer to unity 
were the continuum levels for their spectra drawn higher by 5% - 10%, a possibility raised in §5. 
As matters stand, errors in continuous opacity are limited to <10% nonetheless. 

The apparent V magnitude may also be compared to that expected from the model. From the 
V magnitudes tabulated with the model colors, we calculated the difference in the stellar model V 
with respect to the solar model, added this and 51og(L>/i?) to the observed solar V = —26.7, then 
subtracted this from the observed stellar V magnitude. Again, no reddening or other corrections 
were applied. The resulting differences are shown in column 8 of Table 1. The agreement is 
reasonable: the theoretical magnitudes average 0.09 it 0.06 mag brighter than those observed, with 
a ± 0.18 mag individual standard deviation. 



7. Summary and Discussion 

In brief, we have calculated high-resolution mid-UV spectra from first principles, after modify- 
ing line parameters and assigning line identifications for missing lines based on changes seen among 
spectra of stars spanning a range of well-determined values of Tgfj, log 5, and [Fe/H]. The resulting 
calculations provide a good fit for turnoff stars of all metallicities redward of 2900 A . Previous 
discrepancies in fitting the solar pseudo continuum (noted in §1) are attributed not to missing con- 
tinuous opacity, which is good to <10% near 3100 A, but rather to missing or miscalculated line 
opacity. 
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To improve the accuracy of the mid-UV gf-values, damping constants, and hne identifications 
derived here, we have determined temperatures, gravities, microturbulence, and abundances for 
eight stars by calculating fits to optical as well as UV spectra. For stars more metal-poor than 
one-thirtieth solar, the fit to the mid-UV spectrum was satisfactory using laboratory lines alone, 
and can by itself determine Tefj . For the more metal-rich stars, lines missing from the laboratory 
list compromise the goodness of fit of the mid-UV spectrum. A 100 K increase in T^g is seen to be 
difficult to detect from mid-UV spectra alone, as long as [Fe/H] and log g are increased by 0.1 dex 
to compensate. This Tefj-[Fe/H] degeneracy parallels the age-metallicity degeneracy in elliptical 
galaxies (Worthey 1994, 1999). As with galaxies, it is reduced here by considering the Balmer 
lines, with the advantage that stellar models need not include Balmer-line emission nor composite 
populations. Here Tgff was set for stars with [Fe/H] > —1.5 by matching the Ha line profile and 
the strengths of other high-excitation features. 

For all the metal-poor stars, the excellent match between calculations and observations of both 
optical and mid-UV spectra, and between model and observed colors, indicates an uncertainty of 
±50 K in an individual Tefj determination, and a similar uncertainty in the Tefj scale. Moreover, 
the stellar V magnitudes generally agree to 0.1 mag, and the mid-UV normalization constants to 
10%, with those expected from the stellar angular diameter found from the observed parallax and 
the model log g, assuming reasonable masses. 

In this process, it was found that these optical diagnostics and the slope of the mid-UV con- 
tinuum could only be fit simultaneously when using models of Castelli et al. (1997) in which 
convective overshoot has been turned off. Our work thus agrees with that of Castelli et al. 
(1997), who found that their models better reproduce the Balmer profiles in the Sun than do 
those of Kurucz (1995). Downloading the Castelli et al. (1997) models from the Kurucz web site 
at http://cfaku5.harvard.edu is recommended for all work where the results depend upon the tem- 
perature structure in deep continuum-forming regions. 

Residual discrepancies in fitting the profiles of very strong atomic lines, and the strengths of 
strong versus weak solar OH lines in the 3000 - 3100 A region, are reduced here by mimicking a 
chromosphere towards the surface of every model with [Fe/H] > — 1. As described at the end of 
§4, the model temperature was raised in the region just above the depths where optical continuum 
and weak lines are formed. While this procedure was necessary to match strong-line cores and OH 
line strengths over a very broad range of metallicity with a single set of gf-values, it is definitely 
an oversimplification. Energy and stability considerations remain to be examined. Fits might be 
improved with a different choice of temperature enhancement. Indeed, a two-stream model might 
be found necessary. In the meantime, line cores also will be subject to uncertainties due to the 
presence of chromospheric emission. Peterson &: Schrijver (1997) demonstrate that Mgii emission is 
always present in solar-temperature stars, with a strength that increases with increasing metallicity. 
Emission also increases with increasing stellar activity (Ayers et al. 1995). Since activity in turn is 
associated with rapid rotation, stronger Mgii emission appears commonly among young stars and 
close binaries. Consequently, no single quiescent chromosphere appears capable of modeling the 
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very core of the profile of a strong line in all turnoff stars. 

The reliability of our T^fi values and the good match achieved for OH lines enable a second look 
at the oxygen abundances of metal-poor stars, which are still uncertain as noted in §1. Our values 
for Tcfj are intermediate, suggesting that T^s is not the sole cause of the oxygen discrepancies. 
We defer a more detailed discussion to a future paper by Peterson, Dorman, & Rood (2001), who 
compare oxygen abundances derived from the OH lines with those from the permitted and forbidden 
lines of atomic O i in these stars. 

The manual addition of lines to the model line list, to match otherwise missing features seen 
in the observed spectrum, is clearly necessary for even an approximate representation of the mid- 
UV spectrum of solar-temperature stars of one-tenth solar metallicity and higher. Nonetheless 
the procedure has its important limitations. One is that as overall line strength increases, the 
proportion of missing lines goes up dramatically; disentangling them requires a broader range of 
spectra of high quality, as discussed below. Another is that the excitation or ionization of the 
transition could be misjudged; the ramifications are discussed in the next paragraph. A third is 
that many missing lines may not be due to iron as assumed, but to other elements. As discussed 
in §4, this is potentially serious redward of the bound-free absorption edge of Mgl near 2512 A. 
Elsewhere the effect should be minimal, because although other iron-peak elements contribute 
significant line absorption throughout the mid-UV, the relative abundances of all other iron-peak 
elements with respect to iron change little from one late F or early G main-sequence star to another 
in both halo and disk populations in our Galaxy. To be sure, the list generated by this procedure 
cannot be expected to reproduce the spectrum of a peculiar A or F star, for example a Cr-enhanced 
star, in which the iron-peak ratio has been dramatically altered (Jaschek &: Jaschek 1990). 

More generally, the list is not suitable for modeling stars whose parameters lie well outside the 
range of those analyzed here. In such cases, any difference becomes important between the assigned 
and the actual excitation (and ionization) potential of the transition responsible for a line, leading 
to a significant misrepresentation of its strength. Consequently, the present list cannot be expected 
to reproduce the spectra of stars hotter than early F nor cooler than early G, nor giants, as well as 
stars of metallicity substantially higher than solar for which lines are still missing. The list could 
be made suitable for those stellar types by extending the present analysis to include such stars, 
provided high-resolution, high S/N mid-UV spectra are available in which mismatches induced by 
errors in the list can be visually identified and corrected. We discuss this further below. 

Conversely, the list as it stands should do as well as indicated by Fig. 1 and Fig. 3 in reproducing 
the mid-UV spectra of any other late F or early G star of solar metallicity or less, provided a suitable 
model atmosphere is chosen. This holds because once the gf-value and damping constant of a line 
transition are set, elemental abundance and the Boltzmann and Saha equilibria at each point in the 
atmosphere fully determine line strength in quiescent stars such as these. Dynamical effects can 
be ignored, and non-equilibrium effects (if any) are slowly-varying functions of spectral type, hence 
absorbed into the gf-value determinations. The only exception is the chromospheric emission noted 
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above. Were the line list improved to the point that a reliable match is achieved to the mid-UV 
spectra of F and G stars of solar metallicity and above, we could expect to use the mid-UV slope to 
at least constrain age and metallicity, regardless of its overall level or the appropriate light-element 
abundance ratio, as these could be incorporated in a grid of calculations covering the entire relevant 
parameter space. 

Spectral indexes offer an alternative means by which mid-UV spectral calculations alone might 
constrain metallicity and Tgff of a star or age of a stellar system. The principle is to construct 
indexes based on spectral regions dominated by features whose behavior breaks the age/metallicity 
degeneracy and elucidates the alpha enhancement, examples of which are the Balmer lines and ratios 
of Fe i/Fe ii and Mg i/Mg ii line strengths. Spectral indexes have indeed been used as diagnostics 
of age and metallicity in ellipticals by several of the groups mentioned in §1. For example, Lotz, 
Ferguson, & Bohlin (2000) have compared mid-UV and optical line indexes they calculated from 
Kurucz flux distributions against mean values observed for the stars considered by Fanelli et al. 
(1990). Unfortunately, poor agreement was seen for all but two mid-UV indexes, the slope of the 
pseudocontinuum between 2600 A and 3100 A, and blended absorption near 2538 A. 

We should be able to improve upon some of those index calculations immediately. As Figs. 1 
and 3 illustrate, the current spectral calculations should be reasonably reliable at all metallicities 
in spectral regions redward of 2900 A, and for those dominated by strong lines blueward of that. 
This includes several indexes already in use, notably the Mg ii and Mg i lines at 2800 A and 2852 A 
and line blends in the 3000-1- A region. We plan to begin such index calculations shortly, and make 
them generally available. 

Blueward of 2900 A , however, our calculations are not yet reliable enough. Fig. 1 shows that 
our calculated spectra of solar-metallicity stars substantially overestimate the true fluxes in those 
wavelength regions not dominated by strong lines. Somewhat ironically, our calculations thus 
suggest that neither of well-behaved Lotz, Ferguson, & Bohlin (2000) indexes is currently reliable 
at near-solar metallicity. 

The problem undoubtedly is lines still missing from the line list. As seen in Fig. 3, missing lines 
often overlap in the solar spectrum to such an extent that they could not be individually discerned, 
and so their species and wavelength assignments became arbitrary. This was true especially in the 
2550 - 2600 A and the 2640 - 2710 A regions, where we did not complete the process. This is the 
reason for the excess flux seen in the calculations at solar metallicity in these regions. 

This could be remedied were spectra obtained of the same quality as that of a Cen A for stars 
of the line strength of HD 184499 that span a range of properties, notably Tcfj . For example, the 
spectra of HD 184499 and a Hyades F star with Tcfj ~6500 K should have about the same overall 
line strength: the Hyades star's higher abundance, [Fe/H] ~-|-0.15, is offset by its higher Tgff at the 
rate of 0.1 dex per 100 K, as Fig. 2 illustrates. Because line crowding is much lower than in a Cen A 
, comparing high-resolution mid-UV spectra of HD 184499 and a sharp-lined Hyades F star could 
show where the strongest lines are that are missing at solar metallicity and suggest their identities 
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from the changes observed with Tpg . The a Cen A and solar spectra could then be used to fill in 
the rest. Analyzing an additional star closer in Tcfj to HD 184499 but with [a/Fe] = would show 
explicitly which transitions are due to magnesium in the 2550 A region, and whether OH has been 
properly characterized throughout the 2900 - 3100 A domain. 

At that point, the full arsenal of indexes representing the strengths of mid-UV features could 
be calculated from the spectra, and their behavior with Tefj , log (/, [Fe/H], and [a/Fe] assessed 
explicitly. This is a major advantage of calculations from first principles over empirical libraries, 
which are confined to the spectra of stars available nearby. It is a major goal of this project. 
Indexes should result with significant diagnostic power for the integrated-light spectrum of an old 
stellar system. 
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Fig. 1.— 

Plots are shown comparing observed (heavy hne) and calculated (light line) spectra for the 
nine stars with fluxed mid-UV spectra included in this study. To the right, the name of each star 
in bold is accompanied by the model parameters used for the corresponding calculation. At the 
left is the mid-UV normalization constant for each comparison. 

Fig. 2.— 

Plots as in Fig. 1 are shown comparing the fits for two stars with small changes in model 
parameters. 

Fig. 3.— 

The same as Fig. 1, but on a larger scale in 3 A spectral regions. At the bottom is the Sun: two 
scans of the center of the solar disk are superimposed. At the very top, the strongest lines in the 
solar spectrum calculation are identified as described in the text. The same mid-UV normalization 
constants as in Fig. 1 are adopted for the top nine stars, a) 2627.5 - 2630.5 A. b) 2642.5 - 2645.5 A 
. c) 2900 - - 2903 A. d) 3082.5 - 3085.5 A. e) 3100 - 3103 A. 

Fig. 4.— 

Comparisons as in Fig. 1 are shown for optical spectral regions containing magnesium or 
Balmer lines in eight of the stars plus the Sun. Each observed spectrum (heavy line) has been 
normalized to unity as described in the text. Line identifications as described for Fig. 3 are shown 
for the strongest lines in the solar calculation, a) Mgll. b) Low-excitation Mgl. c), d) The Mg6 
lines, e) High-excitation Mgl. /) Ha. In this panel only, the third plot from the top represents a 
non-optimal choice of model parameters. It illustrates the degradation of the fit to Ha when a T^fi 
is chosen 100 K cooler than that of the plot immediately below. 
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